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We present a general theoretical method to generate maximally entangled mixed states of a 
pair of photons initially prepared in the singlet polarization state. This method requires only 
local operations upon a single photon of the pair and exploits spatial degrees of freedom to induce 
decoherence. We report also experimental confirmation of these theoretical results. 
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I. INTRODUCTION we draw our conclusions in Sec. IV. 



Entanglement is perhaps the most puzzling feature 
of quantum mechanics and in the last two decades it 
became the key resource in quantum information pro- 
cessing Entangled qubits prepared in pure, maxi- 
mally entangled states are required by many quantum- 
information processes. However, in a mundane world, 
a pure maximally entangled state is an idealization as, 
e.g., a plane wave in classical optics. In fact, interac- 
tion of qubits with the environment leads to decoher- 
ence that may cause a pure entangled state to become 
less pure (mixed) and less entangled. Thus, any realis- 
tic quantum-communication/computation protocol must 
cope with entangled mixed states and it is desirable to 
attain the maximum amount of entanglement for a given 
degree of mixedness. States that fulfill this condition are 
called maximally entangled mixed states (MEMS) and, 
recently, they have been the subject of several papers 
(see, e.g., 0, Q and references therein). In this Article 
we propose a new method to create MEMS from a pair 
of photons initially prepared in the singlet polarization 
state. 

Kwiat and coworkers were the first to achieve 
MEMS using photon pairs from spontaneous parametric 
down conversion (SPDC). They induced decoherence in 
SPDC pairs initially prepared in a pure entangled state 
by coupling polarization and time degrees of freedom of 
the photons. At the same time, a somewhat different 
scheme was used by De Martini and coworkers [3] who 
instead used the spatial degrees of freedom of SPDC pho- 
tons to induce decoherence. However, both the Kwiat 
and the De Martini method require operations on both 
photons of the SPDC pair. On the contrary, our tech- 
nique has the advantage to require only local operations 
upon one of the two photons. 

This Article is structured as follows: In the first part of 
Sec. II we show the relation existing between a one-qubit 
quantum map and a classical-optics setup on the labo- 
ratory bench. In the second part of Sec. II, we exploit 
this knowledge to design a simple linear-optical set-up to 
generate MEMS from a pair of photons via local opera- 
tions and postselection. Then, in Sec. Ill we provide an 
experimental demonstration of our method, using entan- 
gled photons from parametric down-conversion. Finally, 



II. THEORY 

We begin by giving a brief description of the connec- 
tion between classical polarization optics and quantum 
mechanics of qubits, as recently put forward by several 
authors [1, Q . Most textbooks on classical optics intro- 
duce the concept of polarized and unpolarized light with 
the help of the Jones and Stokes-Mueller calculi, respec- 
tively [6]. In these calculi, the description of classical 
polarization of light is formally identical to the quantum 
description of pure and mixed states of two-level systems, 
respectively [7|. Mathematically speaking, there is an 
isomorphism between the quantum density matrix p de- 
scribing a qubit and the classical coherency matrix J 
describing polarization of a beam of light: p ^ J/TiJ. 
J is an Hermitean, positive semidefinite 2x2 matrix, 
as is p. A classical linear optical process (as, e.g., the 
passage of a beam of light through an optical device), 
can be described by a 4 x 4 complex-valued matrix A4 
such that (Jout)y = MijM{J\u)kh where, from now on, 
we adopt the convention that summation over repeated 
Latin indices is understood. Moreover, we assume that 
all Latin indices i, j, fc, Z, to, n, . . . take the values and 1, 
while Greek indices a,(3, . . . take the values 0, 1, 2, 3. In 
polarization optics one usually deals with the real-valued 
Mueller matrix M which is connected to M. via a uni- 
tary transformation A : M = K^MK [9*1 . The matrix M 
is often written as 

where (p, d) G M'^, arc known as the polarizance vec- 
tor and the diattenuation vector (superscript T indi- 
cates transposition) , respectively. Note that d is nonzero 
only for dichroic media, namely media that induce 
polarization-dependent losses (PDL) [6i]. is a 3 x 3 
real- valued matrix. It should be noticed that if we choose 
"^00 ~ 1 (this can be always done since it amounts to 
a trivial polarization-independent renormalization) , the 
Mueller matrix of a non-dichroic optical element (d = 0), 
is formally identical to a non-unital, trace- preserving, 
one-qubit quantum map (also called channel) . If also 
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p = (pure depolarizers and pure retarders p), then M 
becomes identical to a unital, one-qubit channel It is 
not difficult to show that any linear optical device that 
can be represented by M (or A/), can also be described 
by a set of at most four distinct optical elements in par- 
allel as M. — AoTq (g) T*, where the four 2x2 Jones 
matrices Tq,, represent four different non-depolarizing op- 
tical elements and Aq, > [1, ■ From the results 
above it readily follows that the most general operation 
that a linear optical element can perform upon a beam 
of light can be written as Jin — + Jout — Tla ^aTaJinTl- 
Since Aq, > 0, the previous equation is formally identi- 
cal to the Kraus form [ij of a completely positive one- 
qubit quantum map £. Therefore, if a single photon 
encoding a polarization qubit passes through an opti- 
cal device classically described by the Mueller matrix 
M = J2a ^aTa T*, its initial state pin will be trans- 
formed according to Pin Pout oc J^a ^oTaPinTl. 

Now that we have learned how to associate a quantum 
map to a set of at most four optical elements, we can 
apply this knowledge to design a simple optical scheme 
suitable for MEMS production. Suppose to have two 
qubits (encoded in the polarization degrees of freedom 
of two SPDC photons, say A and B), initially prepared 
in the state p : p ^ pij^ki\ij){kl\ = ® 
Superscript R indicates reshuffling [13[ of the indices: 
pff. = Pij^ki- Following Ziman and Buzek [iJ] we as- 
sume that p is transformed under the action of the most 
general local (that is, acting upon a single qubit) linear 
map 8 ®X into the state 
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P£=£®I[p](x^KTc,®IpTl®I. (2) 

By writing explicitly Eq. ^ in the two-qubit basis 

= it is straightforward to obtain {ps)ij,ki 

Ea'^«Pmn,ii(^a)™(^a)fe«- Then, from the definition 
of M it easily follows that {p£)ij.ki oc {Mp^)ik,ji- By 
reshuffling pg, this last result can be written in matrix 
form as oc A4p^ which displays the very simple re- 
lation existing between the classical Mueller matrix A4 
and the quantum state p£. Via a direct calculation, it 
is possible to show that if p represents two qubits in the 
sing let state ps = j{I I — ax (Tx — (Jy Uy — az cr^) 
then the proportionality symbol in the last equa- 
tion above can be substituted with the equality symbol: 
P£ = Mpf- Note that this pleasant property is true 
only for the singlet state. However, if the initial state p 
is different from the singlet one, then M must be simply 
renormalized by imposing Tr(A^p^) = 1 . 

Now, suppose that we have an experimental setup pro- 
ducing pairs of SPDC photons in the singlet state ps, and 
we want to transform p^ into the target state pr via a 
local map T(g)I : Ps ^ Pr = i-^Tpf )^- AH we have to 
do is first to invert the latter equation to obtain 

Mr = P^{pf)-\ (3) 
and then to decompose A4t as A4t — J^a^aTa T*. 



Thus, we get the (at most four) Jones matrices Tq repre- 
senting the optical elements necessary to implement the 
desired transformation. This is the main theoretical re- 
sult of this Article. Our technique is very straightforward 
and we shall demonstrate its feasibility later, by apply- 
ing it to design an optical setup devoted to MEMS gen- 
eration. However, at this moment, some caveats are in 
order. To make Air a physically realizable Mueller ma- 
trix, its associated matrix Hq- should be positive semidef- 
inite IJJj]. If this is not the case, then the transformation 
p — > pq- cannot be implemented via local operations. For 
example, it is easy to see that if the initial state is a 
Werner state pw = VPs + ^-^I, (0 < _p < 1) and the 
target state is the singlet pr = Ps, then such operation 
(known as concentration 17]) cannot be physically im- 
plemented by a local setup since Hq- has three degenerate 
negative eigenvalues. Another caveat comes from the no- 
signalling constraint. Since A4t describes a local device 
operating only upon photon A, a second observer watch- 
ing at photon B cannot distinguish the initial state ps 
from the transformed state pq-, that is: p^ = Tr^(ps) = 
Tr^(pr)- This condition requires the one-qubit map T 
to be trace-preserving: J2a ^aT^Ta = I. From Eq. ([1]), 
a straightforward calculation shows that such condition 
cannot be fulfilled if d ^ 0, that is if the device imple- 
menting T contains dichroic (or PDL) elements. 

PDL is important in many commonly used optical de- 
vices as polarizers, circulators, isolators, etc., p. Within 
the framework of quantum information theory, all these 
physical devices may be represented by "unphysicaF one- 
qubit maps T that violate the no-signalling condition. 
This apparent paradox disappears if one allows causal 
classical communications between observers who actually 
measure and reconstruct the target state pr generated 
by the "unphysical" local map T [l^. In fact, in 
coincidence measurements (required to reconstruct pr), 
classical (as opposed to quantum) signalling between the 
two observers is necessary to allow them to compare their 
own experimental results and select from the raw data the 
coincidence counts. In other words, a coincidence mea- 
surement post-selects only those photons that have not 
been absorbed by the PDL element [1]. 

With these caveats in mind, we come to the exper- 
imental validation of our method. We choose to gen- 
erate MEMS I states represented by the density 
matrix pi = p|0+)(0+| + (1 -p)|01)(01|, where |(^+) = 
(|00) + |ll>)/\/2 and (2/3 < p < 1). By varying the pa- 
rameter p, the entanglement and mixedness of the state 
pi change. Here, we use the linear entropy Sl 12011 and 
the tangle T, namely, the concurrence squared [2l|, to 
quantify the degree of mixedness and of entanglement, 
respectively. They are defined as Sl{p) = |[1 — Tr(p^)], 
and T{p) = [max{0, VAq — \Ai — \A2 — n/As}]^, where 
Ao > Ai > A2 > A3 > are the eigenvalues of 
p{(Jy ® ay)p*{ay ® u y) . Aftcr applying Eq. ^ with 
pq- = Pi, a straightforward calculation shows that there 
are only two non-zero terms in the decomposition of A^r , 
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FIG. 1: (color online) Layout of the experimental setup. The 
two-path optical device acts only on photon A. Detectors Da 
and Db perform coincidence measurements. 



namely {Aq 
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2(l-p),Ai {To = 

}. In physical terms, Tq is a polarizer and Ti is 



a 90° polarization rotator. The two eigenvalues {Aq, Ai} 
give the relative intensity in the two arms of the device 
and are physically realized by intensity attenuators. 



III. EXPERIMENT 

Our experimental set-up is shown in Fig.l. Its first 
part (Singlet state preparation) comprises a Krypton- 
ion laser at 413.1 nm that pumps a 1-mm thick (3 — 
BaB204 (BBO) crystal, where polarization-entangled 
photon pairs at wavelength 826.2 nm are created by 
SPDC in degenerate type II phase-matching configura- 
tion |22]. Single- mode fibers (SMF) are used as spatial 
filters to assure that the initial two-photon state is in a 
single transverse mode. Spurious birefringence along the 
fibers is compensated by suitably oriented polarization 
controllers (PC) In addition, total retardation in- 

troduced by the fibers and walk-off effects at the BBO 
crystal are compensated by compensating crystals (CC: 
0.5-mm thick BBO crystals) and half-wave plates (A/2) 
in both photonic paths. In this way the initial two- 
photon state is prepared in the polarization singlet state 
1^^) = {\HV) - \VH))/V2, where H{= 0) and V{= 1) 
are labels for horizontal and vertical polarizations of the 
two photons, respectively. 

In the second part of the experimental set-up (MEMS 
preparation) the two-term decomposition of A4t is physi- 
cally realized by a two-path optical device. A photon en- 
ters such a device through a 50/50 beam splitter (BS) and 
can be either transmitted to path 1 or reflected to path 
2. The two paths defines two independent spatial modes 
of the field. In path 1 a neutral-density filter (Ai) is fol- 
lowed by a linear polarizer (P) oriented horizontally (with 
respect to the BBO crystal basis). When the photon goes 
in this path, the initial singlet is reduced to \HV) with 



probability proportional to the attenuation ratio oi of Ai 
(a — Pout /Pin)- In path 2 a second neutral-density filter 
(A2) is followed by two half wave-plates (A/2) in cascade 
relatively oriented at 45°: they work as a 90° polariza- 
tion rotator. When the photon goes in path 2, the singlet 
undergoes a local rotation with probability proportional 
to the attenuation ratio 02 of A2. 

The third and last part of the experimental set-up 
(Tomographic analysis), consists of two tomographic an- 
alyzers (one per photon), each made of a quarter- wave 
plate (A/4) followed by a linear polarizer (P). Such an- 
alyzers permit a tomographically complete reconstruc- 
tion, via a maximum-likelihood technique [2^ . of the 
two-photon state. Additionally, interference filters (IF) 
in front of each detector (AA = 5 nm) provide for band- 
width selection. It should be noticed that detector Da 
does not distinguish which path (either 1 or 2) a pho- 
ton comes from, thus photon A is detected in a mode- 
insensitive way: This is the simple mechanism we use 
to induce decoherence. In the actual setup, a lens (not 
shown in Fig. 1) placed in front of detector Da focusses 
both paths 1 and 2 upon the sensitive area of the detector 
which becomes thus unable to distinguish between pho- 
tons coming from either path 1 or 2 ("mode-insensitive 
detection" ) . 
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FIG. 2: Experimental data and theoretical prediction (contin- 
uous line) in the linear entropy-tangle plane. The gray region 
represents unphysical states and it is bounded from below by 
MEMS (dashed curve). The lower dotted-dashed curve repre- 
sents Werner states. The horizontal (dotted) line at T = 4/9 
separates MEMS I (above), from MEMS II (below). Stars de- 
note MEMS I states that have the same linear entropy as 
the measured states p™'' (i.e., the experimental points above 
the line T = 4/9). All measured data follow very well the 
theoretical curve. 

Experimental results are shown in Fig. 2 together with 
theoretical predictions in the linear entropy-tangle plane. 
The agreement between theoretical predictions and mea- 
sured data is very good. The experimentally prepared 

initial singlet state has a fidelity ^25|] F{ps,pf^^) — 
2 



cxp 



^) '-^ 97% with the theoretical sin- 



glet state Ps- The continuous curve is calculated from 
the matrix Pc ■ Pc = M-TPf^^i and varying p. It rep- 
resents our theoretical prediction for the given initially 
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prepared state p™^. If it were possible to achieve ex- 
actly pf^P = Ps , then such curve would coincide with the 
MEMS curve above the horizontal (dotted) line T = 4/9. 
Experimental points with T > 4/9 (Pi^^) are obtained by 
varying the neutral-density filters Ai,A2 in such a way 
that 02 > fli; while points with T < 4/9 are achieved 
for 02 < fli. Note that the latter points do not repre- 
sent MEMSs, but different mixed entangled states whose 
density matrix is still given by pi but with the parame- 
ter p now varying as < p < 2/3. The average fidelity 
between the measured states p^^^ an the "target" states 
Pj,, is given by F{p^, p'^'^) ~ 80%. The main reason for 
its deviation from < 100%, is due to spurious, uncon- 
trolled birefringence in the BS and the prism composing 
the set-up. To verify this, first we calculated the fidelity 
between the states Pc{p) (obtained by applying the the- 
oretically determined map T (8) X to the experimentally 
prepared initial singlet state pf^^), with the theoretical 
MEMS pi{p). We have found F[pi{p), pc{p)] > 97% for 
all 2/3 < p < 1; thus the value of ~ 80% cannot be 
ascribed to the imperfect initial singlet preparation. Sec- 
ond, we explicitly measured the Mueller matrices for both 
the BS and the prism (matrices that would be equal to the 
identity for ideal non-birefringent elements) and we ac- 
tually found spurious birefringence. From such measured 
matrices it was possible to determine the unwanted local 
unitary operation induced by these optical elements [1^ . 
It is important to notice that such operation does not 
change the position of our experimental points in the lin- 
ear entropy-tangle plane. Now, if one applies this unitary 
operation to our raw data and calculates once again the 
average fidelity, the result would be _F ~ 91%. However, 
since this "compensation" of the spurious birefringence 
is performed upon the measured data and not directly 
on the physical setup, we felt that it was more fair to 
present the uncorrected data. 

IV. DISCUSSION AND CONCLUSIONS 

In conclusion, we have theoretically proposed and ex- 
perimentally tested a new, simple method to create 



MEMS I states of photons. This method can be easily 
generalized to generate MEMS II states, as well. How- 
ever, this task would require a slightly different experi- 
mental setup with a three-path linear optical device act- 
ing only upon photon A (26j . In particular, we have 
shown that it is possible to create a MEMS from a SPDC 
photon pair, by acting on just a single photon of the 
pair. This task could appear, at first sight, impossible 
since it was recently demonstrated [l^l that even the 
most general local operation cannot generate MEMS be- 
cause this would violate relativistic causality. However, 
as we discussed in the text, our results do not contra- 
dict Ref. [l3] since we obtained them via postselection 
operated by coincidence measurements. The latter are 
possible only when causal classical communication be- 
tween detectors is permitted. Still, the connection be- 
tween relativistic causality, dichroic (or, PDL) devices 
and postselection, is far from being trivial. For exam- 
ple, suppose that a two-photon state is produced by an 
optical setup containing local PDL elements, and that 
we tomographically reconstruct it after coincidence mea- 
surements. Such a reconstructed state will correctly de- 
scribe the result of any other measurement involving co- 
incidence measurements (as, e.g.. Bell measurements), 
but it will fail when describing the result of any single- 
photon measurement. We stress that this limitation is 
not inherent to our scheme, but it is shared by all optical 
set-ups containing PDL elements. 



Acknowledgments 



We acknowledge Vladimir Buzek for useful correspon- 
dence. We also thank Fabio Antonio Bovino for helpful 
suggestions. This project is supported by FOM. 



[1] M. A. Nielsen and L L. Chuang, Quantum Computation 

and Quantum Information (Cambridge University Press, 

Cambridge, UK, 2002), reprinted first ed. 
[2] N. A. Peters, J. B. Altepeter, D. A. Branning, E. R. 

Jeffrey, T.-C. Wei, and P. G. Kwiat, Phys. Rev. Lett. 

92, 133601 (2004). 
[3] M. Barbieri, F. De Martini, G. Di Nepi, and P. Mataloni, 

Phys. Rev. Lett. 92, 177901 (2004). 
[4] N. Brunner, A. Acm, D. Collins, N. Gisin, and 

V. Scarani, Phys. Rev. Lett. 91, 180402 (2003). 
[5] A. Aiello, G. Puentes, D. Voigt, and J. P. Woerdman, 

Opt. Lett. 31, 817 (2006). 
[6] J. N. Damask, Polarization Optics in Telecommunica- 



tions (Springer, New York, 2005). 
[7] D. L. Falkoff, and J. E. McDonald, J. Opt. Soc. Am. 41, 

862 (1951); U. Fano, Phys. Rev. 93, 121 (1954). 
[8] M. Born and E. Wolf, Principles of Optics (Cambridge 

University Press, 1999), seventh ed. 
[9] A. Aiello, and J. P. Woerdman, math-ph/0412061 (2004). 
[10] S.-Y. Lu and R. A. Chipman, J. Opt. Soc. Am. A 13, 

1106 (1996). 

[11] M. B. Ruskai, S. Szarek, and E. Werner, Linear Algebr. 

Appl. 347, 159 (2002). 
[12] D. G. M. Anderson and R. Barakat, J. Opt. Soc. Am. A 

11, 2305 (1994). 
[13] I. Bengtsson and K. Zyczkowski, Geometry of Quan- 



5 



turn States. An Introduction to Quantum Entanglement 
(Cambridge University Press, 2006), 

[14] M. Ziman and V. Buzek, Phys. Rev. A 72, 052325 (2005). 

[15] In this Article the symbols ax,(Jy,Oz denote the three 
Pauli matrices as defined in: J. J. Sakurai, Modern Quan- 
tum Mechanics (Benjamin/Cummings, Menlo Park, Cal- 
ifornia, USA, 1985). 

[16] The Hermitean, positive semidefinite matrix H was first 
introduced in polarization optics by R. Simon, Optics 
Comm. 42, 293 (1982), and subsequently developed in 
Ref. [l^]. It represents the polarization-optics analogous 
of the dynamical (or, Choi) matrix describing a one-qubit 
quantum process [E. C. G. Sudarshan, P. M. Mathews, 
and J. Rau, Phys. Rev. 121, 920 (1961); M.-D. Choi, 
Linear Algebra Appl. 10, 285 (1975)]. 

[17] R. T. Thew and W. J. Munro, Phys. Rev. A 63, 
030302(R) (2001). 



[18] A. Aiello, G. Puentes, D. Voigt, and J. P. Woerdman, 

quanth-ph/0603182 (2006). 
[19] For sake of simplicity, we do not consider here type II 

MEMS. See, e.g., Ref. 3 and references therein. 
[20] S. Rose and V. Vedral, Phys. Rev. A 61, 040101(R) 

(2001). 

[21] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[22] P. G. Kwiat, K. Mattle, H. Weinfurter, A. Zeilinger, A. V. 

Sergienko, and Y. Shih, Phys. Rev. Lett. 75, 4337 (1995). 
[23] G. Puentes, D. Voigt, A. Aiello, and J. P. Woerdman, 

Opt. Lett. 31, 2057 (2006). 
[24] D. F. V. James, P. G. Kwiat, W. J. Munro, and A. G. 

White, Phys. Rev. A 64, 052312 (2001). 
[25] R. Jozsa, J. Mod. Opt. 41, 2315 (1994). 
[26] A. Aiello, G. Puentes, J. P. Woerdman, 

|quanth-ph/0611179| (2006) . 



